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A theoretical description of exclusive nonleptonic decays of charmed hadrons based 
on general principles is not yet possible. Even if the short distance effects due to hard 
gluon exchange can be resummed and the effective hamiltonian has been constructed at 
next-to-leading order , the evaluation of its matrix elements requires nonperturbative 
techniques. A classic analysis based on QCD sum rules has been presented in three papers 
by Blok and Shifman |^], but only the general trends were reproduced: the agreement with 
present data is poor at a quantitative level. Waiting for future progress in lattice QCD 
calculations one has to rely on approximate methods and on models. 

We recently presented one such model, based on the factorized approximation, 
with annihilation terms and rescattering effects due to resonances coupled to the final 
states, that has been rather successful in describing the bulk of the experimental data. Its 
main shortcoming was, in our opinion, the large flavour SU(3) breaking in the axial charges 
forced by the fitting of the data on decay rates to final states with one pseudoscalar and 
one vector meson [PV). 

In this letter we modify the previous approach, inserting rescattering corrections 
also in nonresonant channels. Moreover, we include the decays to final states containing 
one of the lowest mass scalar mesons {S), /o(980) and ao(980), that are connected through 
rescattering effects to the previosly considered PV final states. In this way, we are able 
to obtain a much better fit of the experimental data, while keeping the SU(3) breaking in 
the axial charges at a smaller and more acceptable level. 

The scattering phase shifts in nonresonant channels were neglected in . For decays 
to PP final states this is essentially correct, given that only one nonresonating phase, 
corresponding to the 27 representation, is involved and of course the final rates only depend 
on the phase shift differences. In the case of PV final states, on the other hand, many 
different SU(3) representations are present: to minimize the number of parameters, we 
only include a nonzero phase shift for the 27 (besides the resonant 8^^) and keep the 
others to zero. The 27 phase shift is most welcome, especially to obtain a better fit for 
— > PV Cabibbo-allowed decays. We admit two different values for the phase shift at 
the different energies, corresponding to the masses of D and of Dg: the fitted values are 
527(m£)) = 47.4° and 527(iiidJ = 59°, reasonably similar to each other, although maybe 
larger than expected. 

The nature of the scalar resonances, /o(980) and ao(980), has been discussed for 
quite a long time. They do not look like the members of a normal nonet, in that the fo 
is strongly coupled to KK, and could be for this reason identified with an ss state, but it 
is degenerate in mass with the isovector gq. Moreover, the strange scalar states lie quite a 
bit higher. For these reasons, it has been suggested that the fo and ao are essentially KK 
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molecules and are made therefore of two quarks and two antiquarks, i.e. an ss pair plus 
a light qq pair 0. In charmed meson decays, both /ott"'" and D° fo^s have 

been observed experimentally , , ||] . In the factorized approximation, one would have 
for the decay amplitudes prior to rescattering corrections: 



AUDt ^ /o7r+) = 



G p 



V2 



G p 



V2 



V2 

t^^dC/^s (Ci+eC72)/K </o|a'^(^^)Ml^°>■ 



G'F 



In (|I|) Ci are the Wilson coefficients in the effective hamiltonian, ^ is the color screening 
parameter (that should be equal to if the factorization approach were exact), the 

axial currents are denoted by (A^, = q' '^^'^^q and we neglected possible annihilation 
contributions. The observation of the decay —>■ /ott^ would imply the ss nature for 
/o, while /o-^° points to a nonstrange composition. 

Following the suggestion of we consider /o and ao as cryptoexotic two-quark plus 
two-antiquark states and attribute them to (incomplete) 8 and 1 SU(3) representations 
Wo) e |8), l/o) G xfl\8) + -y/l 11). We then define 



{fo\d''iAl),\Dt)={Mh^-Ml] _ 

(2) 



(fQ\d^(A^,)JD^) ={Ml - Mj )— — 



The axial charge as is a parameter to be fitted. The result is as — 0.39, smaller - as 
expected - than the corresponding axial charges for D transitions to vector mesons. 

Our model also predicts charmed meson decays to states including the ao(980) meson 
and Cabibbo suppressed decays to PS, not yet observed. The amplitudes in factorized 
approximation are easily obtained, and the relevant form factors are all expressed in terms 
of the parameter as'- as an example, 

ialW{A''.)JD+) = -(Ml- Ml )^ — — . (2') 
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We describe now in more detail the procedure followed to include final state in- 
teractions. Defining as B the decay amplitude including the phase space factor, i.e. 
i3w = v^p / (87rm|,) where p is the momentum of the final particles in the D rest 
frame, we have for D —>■ PV (PS) decays 

BiD ^ Vh Pk) =B^{D^ Vh Pk)+Chk[eMi8s) - l\A^T + c^/.fc[exp(z527) - IjAfT , 
B{D ^ Sh Pk) = B^iD ^ Sh Pk)+xps dhkieMiSs) - 1] A| + (3) 

+yps rf/^fc[exp(^527) - 1] Af , 



where 



8 _^h'k' ^h'k' B^{D Vh' Pk') + xps J2h"k" ^h"k" B^{D Sh" Pk' 



— 



.27 _Ylih'k' ^h'k' B^{D — i> Vh' Pk') + yPS J2h"k" dh"k" B^{D — > Sh" Pk") 

— 

Y.h'k' \dh'k'\'^ + Y.h"k" M/i'"-" 1^ 



(4) 



'k' 



In (H) and (H) Chk (dhk) are the PV couplings to 8^ (27), multiplied by a {p/Mp)i 
kinematical factor. This p dependence must be present in the B amplitudes and, as in 
, we include it in the coefficients in order to automatically decouple the channels below 
threshold. The PS couplings to 8^^ (27), multiplied by their kinematical factor (in this 
case (p/Mp)5), are denoted Chk (dhk)- 

We note that the phase shift is determined by the parameters of the resonance 
P appropriate to the decay channel considered (P = i4r(1830) or 7r(1770)), as follows 

r(P) 

sinSs exp(ids) = — ■ (5) 

2{M~-Mo)-iT{P) 

In the isoscalar case, is a free parameter instead. 

The parameters xps and yps are connected with the mixing between PV and PS 
channels. The representations 8^ (for PV) and 8^ (for PS) have the same parity and 
charge conjugation and may therefore naturally mix, xps 7^ 0. The two 27 representations 
have opposite charge conjugation; the zero hypercharge sectors cannot mix if isospin is a 
good symmetry, while the y=±l terms may be mixed with opposite mixing angles yps 
(this is an SU(3) violating effect: SU(3) symmetry requires equal mixing angles for any Y 
value). We required in the fit \yps\ < l^^psl- 

We have to face the problem of enforcing orthogonality between the resonant 8 and 
the non-resonant 27 channels. These would be orthogonal in the SU(3) symmetric limit. 



For the couplings of the singlet parts of /o, and 77 we adopt nonet symmetry. 
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but they are not. For the PS channels, the scalar multiplet is incomplete and therefore the 
orthogonality is badly broken. Even for the PV channels the cancellations that would give 
orthogonality do not actually take place, since we included in the rescattering coefficients 
the kinematical factors. For Cabibbo allowed decays one has J2h"k" Ch"k"dh"k" — 0.16 
^-J^d Y^h'k' Ch'k'dh'k' - 0.008. 

The difficulty may be nicely overcome taking advantage of the mixing between PV 
and PS final states. The orthogonality requirement 

''^Ch'k'dh'k' + xps yps Ch"k"dh"k" = (6) 
h'k' h"k" 

establishes a relation between xps and yps, so that only one of them remains as a free 
parameter. The best fit values are xps — 0.25 and 

{^0.20, for D Cabibbo allowed (doubly-forbidden) decays; 
0.00, for Dg Cabibbo allowed and D first-forbidden decays; 
+0.19, for Dg first-forbidden decays. 

We briefly recall now the aspects of that are not modified in the present approach. 
For the evaluation of the weak decay amplitudes we use the factorization approximation 
and a pole model for the form factors, as in eqs. (p, (^). The weak vector charges are 
assumed SU(3) symmetric: their value, 0.79, is taken from the experimental results for 
D —>■ Keu. For the axial charges we allow some SU(3) breaking, and let them vary in the 
range 0.8 ^ 0.9 independently. The decays to final states including i] or rj mesons have 
been treated following the approach of D'yakonov and Eides [^: the rj-rj mixing angle 
is therefore fixed to -10°. For the decays to PP and PV channels we also consider the 
contribution from annihilation (or VK-exchange) diagrams: the relevant matrix elements 
of the divergences of weak currents are given in terms of two parameters to be fitted, Wpp 
and Wpv, with 



< K-n+\d''{Vf)^\0 >= i (m, - rrid) -j^ Wpp, 

Id 

< K-p+\d^{Ai)^\0 >= - (m, + ma) ^ e* ■ pk Wpv 



(7) 



D 



The final state interactions for the PP channels are dominated by the scalar resonances. 
Only one of them, the strange Kq{1950), has been observed |Ty] in the interesting mass 



region. In ||^ we assumed the existence of a nearby isovector resonance oq and we estimated 
its mass from the equispacing formula 

Ml=M],,-Ml + Ml (8) 
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In the fit we allowed the mass, width and branching ratio in the Ktc channel i of the Kq 
resonance to vary within the experimental bounds. From their best fit values (1930 MeV, 
300 MeV and 63.5%, respectively) we get M^^ = 1870 MeV and T^^ = 299.4 MeV. 

In the nonstrange isoscalar case, only relevant for first-forbidden decays, the 
situation is complicated by the possibility of singlet-octet mixing of not yet established 
resonances. The number of parameters (mixing angles, masses and coupling constants) is 
a priori quite large. We imposed the decoupling of the higher mass resonance from the 
nn channel, which together with the requirement of orthogonality reduces the number of 
new parameters to two: the mixing angle i (p and the difference = m^, — m^^ of the 
mass squared (see for details). Using the fitted parameters, the masses and widths of 
the two scalar isoscalar resonances are (M/^, F/^) = (1789, 354) MeV and (Afj', Tj-') = 
(2127, 328) MeV. 

We performed a least square fit with 15 parameters to the 49 data points or ex- 
perimental bounds for the branching ratios. The results are presented in Tables 1 to 4, 
together with predictions for the channels not yet measured. The values of the eleven pa- 
rameters already used in the previous fits are now: ^ = 0.015, ttcu = cicd = 0.9, = 0.8, 
Wpp = -0.269, Wpv = 0.270, M^- = 1930 MeV, F^* = 300 MeV, r = -0.86, = 47.7°, 
A = 1149.4 MeV and = 236.5°. In the axial charges were Gcu = a-cd = 1-0 and 
Qcs = 0.59, while the other parameters are not changed much. We list again the four 
"new" parameter values: S27{mD) = 47.4°, 527{mD^) = 59°, as = 0.390 and xps = 0.249. 
The values of decay constants, quark masses and resonance parameters not explicitely 
mentioned are identical to the values given in [Q. 

The total is 70.3 (of which 6.2 from two Cabibbo doubly-forbidden decays and 
two decays to PS final states, not included in the previous fits). In ref. 0, was 90 for 
45 data points and 11 parameters. A more detailed comparison of the two fits is shown in 
Table 5. We note that the most remarkable improvement occurs for the PV decays: 

it is mainly due to rescattering in the exotic I = \ channel, that is the only rescattering 
effect present in the Cabibbo-allowed decay amplitudes. The worst single point in 
the fit of ref. 0], the branching ratio B(il>+ K n^) (that was Bj/i = 0.64 % versus 

^ Actually, the parameter to be fitted is the ratio r = gsis/gsss, where gsis is the SU(3) 
invariant coupling of the octet of scalar resonances to a singlet and an octet of pseudoscalar 
mesons and ^888 is the coupling to two pseudoscalar octets Q. Nonet symmetry corresponds to 
r = 1. The branching ratio is a quadratic function of r. 

Denoting by |/o) the lower mass state, we define |/o) = sincp \fs) + cos4> |/i) , |/o) = 
-cos^ l/s) + sincf) l/i). 
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Bexp = 2.2 ± 0.4 %), is now fitted quite well, Bth = 2.47 %. This overcompensates the 
slightly worse fit for the decay Ksp^: ^th = 5.60 % now (5.28 % in 0) versus 

Bexp = 3.3±1.25 %. The greater freedom provided by the presence of the new parameters 
^27 allows the reduction of the SU(3) breaking in the axial constants acu = cicd and Ocs, 
that we imposed not to differ by more than 0.1 in this work. It also allows an apparently 
minor change in the annihilation parameters and in the parameter ^, which now happens 
to be small and positive: this has the effect of improving considerably the success of the fit 
also for the decay D+ Ksn+: Bth = 1-35 % (it was 1.08 %) versus B^xp = 1.37±0.15 %. 
A considerable improvement also occurs for the Cabibbo forbidden decay — » K^K : 
Bth = 0.38 % (it was 0.25 %) versus B^xp = 0.51 ± 0.10 %. 

Concerning the decay rates of Df and D^, the quality of the present fit is comparable 
to the fit in ref. @. In particular, for ^ K rj and —>■ p'^rj the results are 
still unsatisfactory (more than three standard deviations lower than the data points). 
Neither annihilation contributions, nor final state interactions were present for channels 
with positive G-parity and / = 1, like , in 0. In this fit the exotic rescattering affects 
these channels, giving for instance a nonzero branching ratio for the decay un^; 
however, it only slightly lowers (going in the wrong direction) the theoretical prediction 
for Df p~^r] . It might be possible to attribute the discrepancy @ to an annihilation 



contribution, not taken into account here, through the glue components in rj and t) [12 . 

Two out of four data points not included in the fit of are very well fitted, but 
the predictions for the other two are not equally satisfactory. The amplitude for the decay 
D'^ — > foKs is colour suppressed and is further decreased by the rescattering effects in 
our model: the theoretical value is therefore smaller than the experimental datum. The 
doubly-forbidden decay —>■ K~^(p can only proceed through annihilation or rescattering: 
also in this case, the theoretical value is considerably lower than experiment. It should be 



noted, however, that recent data from E791 collaboration [|TB[ do not observe a signal in 



this channel and establish an upper bound slightly less than the central value of E691 in 
reported in Table 1. 

As to the predictions for not yet measured decay branching ratios, the largest among 
them refers to the Cabibbo first-forbidden decay K*~^. The decay amplitude 

is colour favoured in this case, and it has a small interfering annihilation contribution 
instead of the larger, although colour suppressed, contribution present in Cabibbo allowed 
decays. The same is true for the process K^K~^. The rescattering effects 



^ The large branching ratio for Df p^rj is difficult to reproduce in many a model, see 
also 11]. 
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decrease the decay rate for (which is in very good agreement with experiment) and 

increase instead the rate for K*^ . The next bigger prediction, for B(il>^ K^p'^), 
deserves a similar comment: it is also increased (~ 20%) by rescattering effects. Among 
Cabibbo doubly-forbidden decays, we predict the largest branching fractions (~ 5 10""^) 
for the decays — > K~^^*^n^. A check for the assumption we made on the scalar particles 
will be the observation of decays with ao(980) production. The largest prediction for not 
yet observed PS decay channels is B{D~^ GqKs) = 0.32 %. 

We will not present here the predictions for CP violating decay asymmetries, that 
depend strongly on the rescattering phases: therefore, they remain similar to those pre- 
viously published i for the PP final states, and differ appreciably in some cases for the 
PV channels. The largest asymmetries (~ —3 10""^) are now predicted in the decays 
D+ p+T] and ^ ujrj : they are entirely due to exotic rescattering, and were there- 
fore zero in The branching ratios of these decays are however small, so that the best 
candidate should be given by the decays p^n^ and D~ P^t^~ , the predicted 

asymmetry being approximately —2 10""^. 

A considerable interest has been recently devoted to the interplay of — mix- 
ing and doubly Cabibbo forbidden amplitudes in the time evolution for decays |1T5[] . 



Particular attention has been given to a term proportional to AM and providing linear 
correction to the exponential decay, present as a consequence of CP violation and/or final 
state interactions, as a possible source of information on "new physics" . A term propor- 
tional to AF is also present. The short distance contributions predicted by the standard 



model are very small for both AM and AF [16]. It was suggested that the mixing may be 



dominated by long distance (hadronic) contributions |]T^ that could result in mixing pa- 
rameters X = AM/F and y = AF/(2 F) as large as 10~^, although this was later criticized 

In our model, we can make an estimate of the long distance contribution to AF 
coming from the two-body states that we included in our fit. This quantity should vanish 
in the SU(3) limit, through an exact cancellation of the contribution of Cabibbo allowed 
and doubly-forbidden transitions with the contribution of once-forbidden decays [0. In 



the presence of SU(3) breaking the cancellation is however not complete. As a consequence, 
our prediction for AF is subject to a large uncertainty; on the other hand, it is to be noted 
that the prediction is independent on the rescattering, provided that, as we impose, the 
sum of the branching ratios remains the same before and after rescattering corrections. 



^ We remark that all the asymmetries reported in Table V of ref. in correspondence to 
D° decays have a wrong sign. The signs for the charged D decay asymmetries are correct. 
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We have {\D^) = CP \D^)) 

= ^ /) ^ /) ^ (i.5 + z0.ooi4)io-^ri5o (9) 

l/> 

In the sum has been approximated including only the contributions of PP (2/3), PV 
(1/3) and PS (~0) final states. Note that the contribution to ri2/r£)0 coming from 
Cabibbo first-forbidden decays alone is 35.2 10"'^, showing that the SU(3) cancellation 
is still rather effective. Although larger than the short-distance prediction, our estimate 
is much smaller than the present experimental bound \y\ = \Ti2\ / T^jo < 0.08. The 
positive sign of the real part of means (if taken seriously) that the shorter lifetime state, 
Dg, decays dominantly into CP-even final states, similarly to the neutral K mesons. 

M.L. acknowledges an interesting dicussion with L. Maiani, S. Malvezzi and D. 
Menasce. 
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f^ 






h 


B,,p(D+ ^ /,) 




KsTT^ 


1.37 ±0.15 


1.35 




0.25 ±0.07 


0.19 







1.70 




0.75 ±0.25 


0.34 


K*\+ 


2.2 ±0.4 


2.47 




< 0.9 


0.73 


Ksp+ 


3.30 ± 1.25 


5.60 


lfK+ 


0.78 ±0.17 


0.81 


KlP+ 





6.30 




< 0.14 


0.13 


at Ks 





0.32 







0.44 


ttQ Kl 





0.24 




< 1.2 


0.013 







0.056 




< 1.5 


0.12 


K+T] 





0.018 




< 0.7 


0.019 


K+ri' 





0.031 




0.67 ±0.08 


0.61 







0.019 







1.71 







0.048 


k*'k+ 


0.51 ±0.10 


0.38 


K*+ri 





0.030 


/o7r+ 


— 


0.028 


K*+r]' 





0.0002 







0.059 




— 


0.030 




— 


0.012 


K+uj 




0.021 


4 V 




0.074 


K+(p 


0.039 ±0.022 


0.0051 


K+fo 




0.0023 












0.0062 



TABLE 1 

Branching ratios for nonleptonic decays. 
[Experimental data and 90% c.l. upper bounds from ref. 
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KsK+ 


1.75 ±0.35 


2.37 


K+n^ 





0.14 


KlK+ 





2.09 


K+r] 





0.28 




1.90 ±0.40 


1.23 


K+ri' 





0.44 




4.7 ± 1.4 


5.39 




< 0.7 


0.40 


p+r] 


10.0 ±2.2 


7.49 







0.044 


p+r] 


12.0 ±3.0 


2.41 


K+p"^ 





0.29 




3.3 ±0.5 


3.96 


K*+ri 





0.18 


KsK*+ 


2.1 ±0.5 


1.87 


K*+r]' 





0.025 


KlK*+ 





2.13 


K+uo 





0.15 




3.5 ±0.4 


4.08 


K+(p 


< 0.25 


0.018 




< 1.7 


0.26 







0.29 


p^n+ 


< 0.28 


0.24 


K^p+ 


— 


1.39 


p+n"^ 


— 


0.24 


foK+ 


— 


0.069 


fon+ 


1.0 ±0.4 


1.06 


a+RO 




0.003 


^0 V 




0.007 


aO K+ 




0.007 


+ ' 




0.002 


K*^K+ 




0.008 



TABLE 2 

Branching ratios for Df nonleptonic decays. 
[Experimental data and 90% c.l. upper bounds from ref. 0] 
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f^ 
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Bexp(I^° ^ fi) 




K-TX+ 


4.01 ±0.14 


4.04 







0.052 


KsTT^ 


1.02 ±0.13 


0.72 


' 
TV r] 





0.16 







0.53 


rjrj 





0.088 


Ks rj 


0.34 ±0.06 


0.42 


r\r\ 





0.18 


KlV 





0.31 




0.088 ± 0.023 


0.110 


Ks T]' 


0.83 ±0.15 


0.78 




0.159 ±0.012 


0.159 


Kl v' 





0.61 


K+K- 


0.454 ±0.029 


0.446 


K*\^ 


3.0 ±0.4 


3.49 


Km"" 


0.11 ±0.04 


0.098 


Ksp^ 


0.55 ±0.09 


0.47 







0.014 







0.33 







0.020 


K*-TT+ 


4.9 ±0.6 


4.85 


' 
p T] 





0.008 


K-p+ 


10.4 ± 1.3 


11.02 


ujrj 





0.20 


K rj 


1.9 ±0.5 


0.37 


urj 





0.0001 


K rj 


< 0.11 


0.004 







0.11 


Ks u 


1.0 ±0.2 


0.88 







0.090 


Kluj 





0.80 




< 0.08 


0.064 


Kscf) 


0.415 ±0.060 


0.40 




< 0.15 


0.062 


Kl(I> 





0.42 


K*+K- 


0.34 ±0.08 


0.43 


foKs 


0.23 ±0.10 


0.037 


K*-K+ 


0.18 ±0.10 


0.30 


foKL 


— 


0.031 




— 


0.69 


a^oKs 




0.109 






0.57 






0.083 






0.12 


4 K- 




0.078 









TABLE 3 

Branching ratios for Cabibbo allowed and first-forbidden decays. 
[Experimental data and 90% c.l. upper bounds from ref. 0] 
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f^ 






/. 











0.0006 


K+n- 


0.031 ±0.014 


0.033 


fo r] 





0.004 







0.0039 




— 


0.011 


K*+n- 


— 


0.035 


«o V 


— 


0.015 


K+p- 


— 


0.025 






0.003 






0.009 


7r+ 




0.070 


K*^r]' 




~ 10-5 












0.004 



TABLE 4 

Branching ratios for Cabibbo first- and doubly-forbidden decays. 
[Experimental data from ref. Q] 



Decays 


7^ data 


(ref. i) 


(This work) 


D+ PP 


5 


9.56 


5.34 


D+ PV 


8 


29.55 


8.46 


D+ ^ PP 


4 


8.79 


7.10 


D+ ^ PV 


8 


15.35 


17.64 


PP 


8 


8.44 


8.43 


PV 


12 


18.35 


17.17 



TABLE 5 

Comparison of our results with the fit of ref. . 
Only Cabibbo-allowed and first-forbidden decays are included. 
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